Gaussian Warp Factor: Towards a Probabilistic Interpretation of Braneworlds 
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We investigate Gaussian warped five-dimensional thick braneworlds. Identification of tlie gravi- 
ton's wave function (squared) in tlie extra-dimension witli a probability distribution function leads 
to a straigfitforward probabilistic interpretation of braneworlds. The extra-coordinate y is regarded 
as a Gaussian-distributed random variable. Hence, all of the field variables and operators which 
depend on y are, also, randomly distributed. Four-dimensional measurable (macroscopic) quantities 
are identified with the corresponding averaged values over the Gaussian distribution. The present 
scenario represents a new phenomenological approach to smooth thick branes which can not be 
obtained through 'smearing out' Randall-Sundrum-like (thin) braneworlds. 
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I. INTRODUCTION 

Braneworld models [H-[3] represent an interesting al- 
ternative to Kaluza-Klein (KK) compactification even if 
the extra-space can be of infinite extension Q • In spite 
of the success of braneworlds models with a 0-width in 
the extra-space, at least at string scale the thickness of 
the brane cannot be neglected. The thin-brane approxi- 
mation is valid as long as the explored energy scales are 
much smaller than the inverse thickness of the brane. In 
for instance, the width of the brane was assumed to 
be of the order of the distance at which the electroweak 
interactions have been probed ~ tti'^]^ {nriEw ~ 10'^ GeV 
is the electro- weak energy scale) . These more realistic al- 
ternatives to 0-width brane configurations are known as 
thick braneworlds A problem with these models is 

what to identify as a four-dimensional observable quan- 
tity. In cosmological applications, for instance, due to 
the finite thickness of the brane, there is some arbitrari- 
ness in the definition of what the effective 4D quantities 
should be. The simplest prescription one can envisage 
is to define the 4D effective quantity associated to a 5D 
quantity as its spatial average over the brane thickness 
[iq | (see also [lil). Needless to say that the above pre- 
scription is legitimate only if we deal with branes with 
well-defined thickness. However, in the case when the 
brane is smoothly spread over the extra-space, there is 
not a unique way to choose the brane thickness, and this 
leads to ambiguities in the computation of 4D measurable 
quantities. 

In the present paper we shall propose an alternative 
way of computing effective quantities within 5D smooth 
thick brane contexts without ambiguities, even if the 
brane thickness can not be chosen in a unique way. Our 



approach will rely on the following assumptions: i) the 
background geometry is a warped one with the warp fac- 
tor being a Gaussian in the extra-coordinate, and ii) the 
extra-coordinate is to be regarded as a continuous ran- 
dom parameter with a Gaussian probability density. The 
resulting picture amounts to a probabilistic representa- 
tion of braneworlds. In this approach the large hierarchy 
between the TeV and Plack scales may be explained if as- 
sume an additional chameleonic interaction between the 
matter degrees of freedom living in the brane and the 5D 
scalar field. 



II. THE MODEL 

In this paper we shall explore thick braneworlds sup- 
ported by a scalar field. The thick branes are described 
by the following 5D action: 



^5 = 



(1) 



where k\ ~ [J^h is the 5D Planck mass), and V{ip) 

is the scalar field's (</?) self-interaction potential. Besides, 
(Vs^/f)^ = g^"^'P,Mf,N , O^ip = g^'^'^Lp-MN- Here cap- 
ital Latin indexes A, B, M, N... — 0,1,2,3,5, while 
small Greek indexes a, f3, /i, i/, ... — 0, 1, 2, 3. The de- 
rived 5D field equations are the Einstein's field equations 
(EFE), and the Klein-Gordon equation (KGE) for the 
scalar field: 



(2) 



respectively, where Gab = Rab — gABRs/^ is the 5D 
Einstein's tensor, and 
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Tab = f,A(p,B - ^5ab(V5(/9)^ - gABV, 



(3) 
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is the scalar field's stress-energy tensor. 

To start with, let us to assume a regular 5D back- 
ground metric which respects 4D Poincare invariance: 



dsg — a'^{y)rif^^dx^''dx'^ + dy^, 



(4) 



where rjajs is the 4D Minkowski metric, a^(y) is the warp 
factor, and, for the moment being, only dependence on 
the extra-coordinate y is being considered. For definite- 
ness we will assume that all of the coordinates x = x^^, 
y, are measured in the same length units, while the warp 
factor is dimensionless. In view of (Hj), the field equa- 
tions ([2]) can be explicitly written in the following way: 



7^ 



2V) , 



6i/2 ^ _A (<^/2 _ 2V) , ip" + AHip' = dV/d^, (5) 

where H = a' /a, and the tilde denotes derivative with 
respect to the extra-coordinate y. Only two of the above 
equations are independent from each other. 

Due to the additional degrees of freedom associated 
with the scalar field (and its self-interaction potential) 
one may choose a suitable functional form of the warp 
factor from the start and, then, one solves the equations 
that involve the scalar field and the self-interaction po- 
tential. In the present case, in order to derive exact so- 
lutions, we make the following (Gaussian) ansatz for the 
warp factor: 



The potential V{(p) in ([7]) makes the given scalar field's 
solution unstable. This has no consequences as long as 
(fi = (p{y) and minimal coupling of matter is considered, 
since the 4D gravity and the evolution of the scalar field 
are decoupled in this case. However, if consider addi- 
tional non-gravitational interaction of matter with tp, the 
instability induced by the above potential can have catas- 
trophic consequences. Nonetheless, if assume ip = </?(x, y) 
from the start, this problem is overcame as a result of the 
chameleon effect (see section WJ^ . 

A few additional comments on our solution. The RS 
solution Q can not be recovered from (|6|7|) in the limit 
A — >• (^ — > oo). Hence, our model does not repre- 
sent a thick brane generalization of the RS model as 
it is the case for a large body of thick brane scenarios 
found in the bibliography [S-Q. It is just a different phe- 
nomenological scenario. Notice that p{y) in ([7]) has not a 
kink-like profile as it is customary in the bibliography on 
smooth thick branes made out of a scalar field (see, for in- 
stance, Ref.jsl). The fact that the solution (|6I7|) supports 
a smooth thick brane configuration will be evident below 
when we study the gravitational content of our model. 
It will be seen that a massless bound gravitational state 
does actually exists which is localized around the origin of 
the extra-coordinate. The smooth thick brane is shaped 
in the extra-space by the wave-function (squared) of this 
gravitational state which can be identified with the 4D 
graviton. 



III. GRAVITON'S WAVE FUNCTION 



a^iy) = cxp(-^^y2). 



(6) 



where the parameter /i is the inverse of the brane thick- 
ness A = As long as a'^{y) = a?{—y), here we shall 
be focusing in 5D solutions which respect mirror symme- 
try: p{y) = ±ip{-y), V{p) = V{~p). 

We can combine equations ([5]) into derived equations 
that can be easier to handle: ZH' = —n^p'"^ f{y) ~ 
± -\/3/i^/k| y -f C, where the " ±" signs account for two 
possible branches of the solution, and C is an integration 
constant which has to be set to zero if one requires mirror 
symmetry to be a symmetry of the solution. Besides, 
the second equation in ([SJ can be rewritten as, k\V = 
kIp>''^/2 - GH"^, i. e., V = 3/^2(1 - 4:^l^y'^)/2Kl Hence, 
we obtain the following solution for the scalar field ip and 
the potential V{ip): 



Here we shall consider linear perturbations of the met- 
ric gl): dsl p^^t = a?{y) [77^^ -t- e /i^j.(x, y)] dx^'dx'' + dy^, 
where ehf^^(:K,y) are small perturbations around the 
Minkowski metric. We shall consider the transverse and 
traceless (TT) gauge, h — rj^^'^hf^i, = rj^^ha^.i, ~ 0, since, 
in this gauge, the linear perturbations of the metric de- 
couple from perturbations of the scalar field 6] . The lin- 
early perturbed EFE-s read: a~^n(^)/i^^-f /iJ^j,-|-4i//(,^^ — 
2{3H' + 8H^)hf,^ = nliip'^ + 2V)hf,^, or, ff make the 
ansatz /i^i,(x,y) = a^'^{y) '^{y) 7r^i/(x), where ^{y) is the 
wave function of the gravitational modes in the extra- 
dimension, then the latter equation can be expressed in 
the form of a Schroedinger equation. 



-Q:^ + UQM{y) 



m a 



(8) 



^(y) - ±yv7^ y, V{^) = - 2^.2 ^2 



(7) 



where Vq = /2k\. This solution is characterized by 
diverging values of the curvature invariants as y — > 00. 
Nevertheless, as it will be discussed in section IIV| this 
does not place any problems to the picture we are about 
to expose. 



where UQM{y) — ^2fi^{l — 2/i^y^) is the Schroedinger- 
like quantum mechanical potential, and accounts for 
the 4D mass of the excited KK gravitational modes: 
□ (^)7rpi,(x) = m^7rpy(x). Equation ([8]) can be explicitly 
written as 



N. 



(9) 
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tinuum of massive states (corrections):^ 



FIG. 1: Probability density profile for the gravitational states. 
The 4D graviton's probability density '^o{y) (solid line), 
and the probability density for the continuum of massive 
states '^miv) (mass parameter m = 0.01- dashed curve, and 
m = 0.013 - dotted line) are shown in the left-hand figure. 
The brane thickness has been arbitrarily set A = 0.5. It is 
seen that the lightest mode of mass m = 0.01 is localized 
farther away from the brane. In the right-hand figure the 
4D graviton's probability Vl'o(y) is shown together with the 
Schroedinger-like potential Uqm{v) times 10~^ (doted curve). 



For the massless mode (properly the 4D graviton): 
2^2(1 - 2^'^y'^)'^Q = 0, which solution is ^"0(2/) = 
sp (— /i^j/^), where, if normalize the wave function 
to unity, )^^^l{y)dv = 1^ N = ^2j^ ii. The 
solution of Eq. ^ for the continuum of KK states is given 
by: ^m{y) = ^/V/2{m/2^lf^o{y) E{y), where 



E{y) = -erf(^?/) erfi(%/2/xy) - 



erfi(e) 



and, erf(a;) = 2 e~^^ d^/^/^:, erfi(a:) = 2 e^^ d^/^/n, 
are the error function, and the imaginary error function 
respectively. 

In the figure [1] we show the profiles of the probability 
densities of the graviton ^o(y) ^'^"^ '^^ the massive states 
4'^(j/). The massless graviton is localized on the brane 
since the analog quantum mechanical potential Uqm is 
well-shaped with a negative minimum inside the brane 
Q. It is seen that the continuum of massive states is 
localized away from the brane which is shaped by the 
massless graviton's wave-function. Besides, the lightest 
of the two massive gravitational modes shown (mass pa- 
rameter m — 0.01) is localized farther away from the 
brane than the heaviest one (m = 0.013), i. e., heavier 
modes are closer to the brane. 



Corrections to the Newtonian Potential 



v{r,y) ~ ^ ^-o + ^i 



-^-^ *i dm. 



where in the second term in the right-hand-side (RHS) 
one has to go to the limit e — 0. After integrating within 
the correction's term we get 

, . K?mim2 ^9, , 37rK?mim2 -,9, ,^9, , 
r r° 

Next we integrate over an effective width of the brane 
T to get a y- independent gravitational potential: 



, , GNmiin2 , 3ttj 1 
VN[r) ^ ( 1 



4/^4 rS J ' 



(10) 



where Gn = ^'^(y) dy, and 

7 - / ^{y)E^{y)dy/ / ■^l{y)dy, 



is a small constant since the massive gravitational states 
are localized away from the brane as seen in the figure [T] 

We have seen that there is a normalizable bound zero- 
mass gravitational state whose wave-function shapes the 
form of the brane in the extra-space. The continuum 
of massive KK modes produces only very small (negli- 
gible) corrections to the Newtonian gravitational poten- 
tial which fall-off very quickly as oc . This is ex- 
pected since the analog quantum mechanical potential 
Uquiy ±c») > (in fact it uncontrollably grows up 
at large |2/|-s). In this case the excited KK-states are 
separated by a gap from the ground state In our 
set-up this may be understood in the following way. The 
massive modes are localized away from the brane: the 
lighter ones are farther away while the heavier KK states 
are closer to the brane. Hence, the lighter states may not 
be excited neither by low-energy processes living in the 
thick brane, nor by self-gravitating interactions with the 
bound graviton, since the probability to find the latter 
far from the brane goes like cx exp (— 2/j,^?/^) and van- 
ishes very quickly. Only heavier KK modes which are 
close enough to the thick brane might by excited by in- 
teractions with the graviton, but this would need of an 
energy of the order of the KK mode's mass, and so are 
effectively decoupled from the thick brane also. 

A crude estimate of the effective mass gap mg// can 
be given if make the following reasonable assumptions: 
i) there is a position in the extra-space where the 
contributions from the massless and excited KK modes 



The j/-depending gravitational potential between two 
point-like sources of masses mi and m2 respectively, 
which are located at ?/, arises from contributions of the 
massless graviton 4'o (standard piece) , and from the con- 



This equation is a slight modification of the one in Ref. 7), to 
include the possibility for the point masses to be located at dif- 
ferent positions within the thick brane. 
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are of the same order ^'o(?/*) — ^'^(j/*), and ii) E'^{yt,) 
is of order unity. Under these assumptions rrif.ff ^ 
2(2/7r)i/4/i - (2/7r)i/42/A (recall that A is the brane 
thickness). Hence, for sharper probability densities the 
effective mass gap is larger. If assume, for instance, 
that the width of the brane is of the order of the dis- 
tance at which electroweak interactions have been probed 
~ m^Jy , then the effective mass gap is of the order of the 
electroweak energy scale ~ tuew ~ 10'^ GeV. 



IV. GAUSSIAN PROBABILITY DISTRIBUTION 

After the Gaussian ansatz (jG)) for the warp factor, the 
next step is to regard the extra-coordinate y as a random 
variable which follows a normal (Gaussian) distribution 
centered at the origin y = with probability density 
function 

f{y) = e ^> / I{v)dy^l, 

\l At: a J-oo 

where is the standard deviation. Since the massless 
gravitational state's wave- function ^o(?/) shapes the pro- 
file of the brane in the extra-dimension, it seems quite 
reasonable to identify the graviton's probability density 
^o(y) with the probability density function for the ran- 
dom variable y: 

f{y) ^ ^^,l(y) = ^le~'^, (11) 

where, for definiteness, we identify 2a with the brane's 
width A — 1/fi. This way one endows the braneworld 
picture with a probabilistic interpretation where all of the 
field variables which depend on y behave randomly and 
are Gaussian-distributed quantities. Besides, the prob- 
ability distribution function ^'□(y) induces a "natural" 
prescription for what to identify as a four-dimensional 
observable within thick brane contexts. Actually, in the 
present case a 4D measurable (effective) quantity asso- 
ciated to a 5D quantity Q(x, y) coincides with its aver- 
aged value over the Gaussian distribution: ((5(x, y)) = 
jr^Q{^,y)^l{y)dy.' 

For illustration let us to compute the averaged 5D Ricci 
curvature: (R^) = Hence, even if i?5 diverges at y- 
infinity as oc — y^, its averaged value {R5) - the one mea- 
sured by a 4D observer - is a finite quantity. This means 
that the 5D gravity, as understood by a 4D observer, is 
asymptotically de Sitter. 



^ Integrating over the infinite extent of the extra-space is legiti- 
mate only for processes which occur at energies below the scale 
at which the KK modes couple effectively to the 4D graviton. 
Otherwise one should integrate over an effective brane thickness. 



A. Macroscopic Quantities and Probabilities 

Under the above prescription for a measurable 4D 
quantity, the thick brane picture acquires a simple proba- 
bilistic interpretation. Given that the extra-coordinate y 
is regarded as a random variable, then, the y-dependent 
gravitational and matter fields are submitted to 'micro- 
scopic' random field variables /ii(x, y), such as the met- 
ric g^iy(x, y) (and derived curvature objects), the stress- 
energy tensor of matter T^™^(x, y), etc. Their corre- 
sponding 'macroscopic' (measurable) quantities Mi(x), 
are computed through Gaussian averaging: 

/oo 
M^(x,y)^'2(y)dy. (12) 
-00 

One can average microscopic operators li, equations, 
etc., and, in particular, the average of the microscopic 
field equations = 0, leads to the corresponding macro- 
scopic equations, i. e., those which represent 4D physical 
laws testable by a 4D observer like us: Ei = (ci) = 0. 

One can also introduce probabilities for measured 
quantities as it follows. The 'probability' V of measuring 
fii within the y-interval [y, y+Sy] around the position y — 
yo, is proportional to {V^^Jy^ oc ^^(x, yo) ^'o(yo) (^y = 

■\/2/7rexp (— 2yQ/A^) /ii(x, yo) Sy/ A. In particular, the 
ratio of the probabilities to measure fj,i around y = ±yo 
and around the origin y = (same interval Sy) is 
i'P^. )y„ I (^mJo = exp (-2yoVA2) ^,(x, yo)//i,(x, 0). 



V. EFFECTIVE 4D GRAVITY 

In standard Kaluza-Klein [l2|, [l^l as in the alterna- 
tive Randall-Sundrum (RS) compactification schemes 
after ensuring that a bounded massless graviton exists 
and can account for recovering of Newtonian gravity, 
one needs also to get a 4D effective description of the 
laws of gravity. This is usually done by straightfor- 
wardly integrating in respect to the extra-coordinate(s) 
under the action integral over the extension of the extra- 
dimensional manifold. In this paper we propose an al- 
ternative scheme. Here a bounded graviton exists by the 
same reason it does in RS compactification [3]: because 
of the warped property of the metric (U). However, un- 
like in standard KK and RS procedures, in our approach 
one does not integrate over y within the action ([l} to get 
to an effective 4D picture. Instead, one regards the 5D 
microscopic field equations ^ as fundamental, and then, 
the macroscopic equations - those which describe the 4D 
laws of gravity - are obtained by computing the Gaus- 
sian average of ([2]) . In this and the following sections we 
shall discuss our approach in more detail. 

Let us consider the metric given by the line element 

dsl = a^{y)g^,{x)dx''dx'' + dy\ (13) 
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with the Gaussian warp factor defined as in where 
in the line element (U) we have replaced the Minkowski 
metric by an arbitrary metric, 77^11/ — >■ ^^^(x). In this 
case the 5D 'microscopic' Einstein's field equations plus 
Klein- Gordon equation ([2]), can be written as: 



G 



3H' + 6H^ + ^{ip'^ + 2V) 



1 

= nlip'ip^f, , a~^aip + if" + AHif' = dV/dip, (14) 

where, as before, H — a' /a, and the quantities Qi are de- 
fined in terms of the 4D metric 5^1,: = 
R = g^^R,.., (Vv?)2 = g^'^if^^,^,,, Dip = The 
second equation in dH]) comes from G^fj, = R^^ = 0. 
This equality forces <p(x) = (^oi i- e., the scalar field can 
not depend on the 4D spacetime point. Without loss of 
generality we can choose (p = ip{y). After this, and re- 
calling that a?{y) is given by the Gaussian function ([SJ, 
the microscopic field equations are: 



Gai 



T 



(m) 



AB 



T 



AB 



, D^ip = dV/dip, 



(15) 



where 




/5g^^ , is the mat- 



ter stress-energy and is the one for the scalar field 



(|3]). In terms of the metric (fT3|) equations (|T5]) are trans- 
formed into the following system of equations: 



2 rp{m) 
^5 J /ij/ 



G tiu 



3iJ' + 6ij2 + -^(/ 



2V) 



R 
2^ 



^2 rpim) 
'^S 55 



6H^ 



2V) 



T, 



(m) 



/j5 



ip^^tp' = 0, ip" + iHip' = dV/dip. 



(16) 



As seen, if one requires that Lp — ip{y) tp^^j, = 0, then 
it has to be assumed that the 5D matter stress-energy 
tensor does not have mixed components: T^™"* = 0. Be- 
sides, since we have assumed (|6l7p is valid, then from 

(HSl) it follows that: G^^ = ^1 T^'"\ 4T^T^ = -R/2a'^. 
These equations lead to the following conditions on the 
5D matter's stress-energy tensor: 



G 

a 

R 

2^ 



3H' + 6H^ 



2V) 



dV_ 

dip 



One can see that the solution (|6l7p amounts to the 
macroscopic 4D vacuum EFE-s G^^ = (i? = 0). Where 

the average of the Einstein's tensor ^G^i,^ = Gfj.iy, since 

these quantities do not depend on y, and, for quantities 
which do not depend on the extra-coordinate (including 
constants) the Gaussian average coincides with the given 
quantity: (P(x)) P(x). 



VI. GRAVITATING MATTER 

One may wonder under which conditions the solu- 
tion dH), d?]) is still valid and, at the same time, 4D 
gravitational equations other than those for vacuum 
are obtained for the background metric (IT3t . In or- 
der to seek for an answer one is to add a matter piece 
cx / d^Xy/lg^lCmiXy ^Xi9ab) to the action ((1]), where 
C,n is the Lagrangian density of the matter degrees of 
freedom x- The microscopic field equations obtained are 



T, 



(m) 



M5 



0, Tk 



(m) 



(17) 



If these conditions are met, then standard 4D microscopic 
EFE-s: G^u — i^l 1)1™^ (x, y), are recovered from Ea.([T5|). 
These are to be regarded as a set of one-parametric par- 
tial differential equations. The continuous parameter y is 
what makes the microscopic T)l™''-s to behave as random 
quantities. 

It is difficult to find a 5D stress-energy tensor derived 
from a Lagrangian which meets the conditions Eq. (jl7p 
at once, notwithstanding, here we assume it is possi- 
ble in principle and defer this subject for future work. 
A trivial example of a stress-energy tensor which meets 
(1171) is supplied by 4D matter in the form of radiation: 



9 ^ til' 



^ 



(rad) 



T, 



{rad) 



55 



0. In this case we can write 

Prad {4:U^u^ + g^y) /3, where prad = Pra<i(x) is 
the energy density of radiation measured by a 4D ob- 
server. Also, since = au^, = u^/a, g^^, — a?g^v, 



then T, 

rfiirad) 



(rad) 



2rfi{rad) 



a^Prad {4:U,tUp + g^„) /3, where 
r^l/""'' is the macroscopic stress-energy tensor of the ra- 
diation. The macroscopic Einstein's equations obtained 
after averaging the microscopic ones are 



Gu 



rp[rad) 



SttG 



N 



rp{rad) 



SttG 



N 



^ Needless to say that the semicolon in this last equation refers to 
covariant derivative in respect to the 4D metric g^jiy. 



bmce Af|i = .Jij2/nl, then, the mass hierarchy can not 
be addressed if minimal coupling of matter to gravity is 
assumed. 
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FIG. 2: The chameleon potential 14// l|22[) is shown for the 
following arbitrarily chosen values of the parameters: Vo = 1, 
/i = 1, /3 = 0.5, and r — 0.1 - thick solid curve, t — 5x 10~^ - 
thin solid line, r = 2.5 x 10^^ - dotted line, r = 10~^ - dashed 
curve. In general r — r(x) is a function of the spacetime point 
so that the 4D dynamics affects the position (and existence) 
of the minimum. 



A. Chameleon Coupling of Matter 

In this subsection, in order to be able to address the 
mass hierarchy issue within our approach, we shall allow 
for a non-minimal (chameleon) coupling of matter to the 
scalar field where, in harmony with effective string theory 
[l3|, the strength of the coupling /3 should be of order 
unity. In this case it is fundamental that the scalar field's 
potential V{(p) be stable, since otherwise, there will be 
a serious problem with stabilization of the y-position of 
the Gaussian distribution function for the matter. 

Here we assume ip = (^(x, y) = ip{y) + (/'(x). Hence, 
since Gg^ = 0, g^'^^'^^^+'^^rg^"^ = 0. We 

shall consider a matter field x which couples to the 
conformal metric in the way a chameleon does 



gAB 



gAB- Hence, the chameleonic matter piece 



action is given by Si^ 



(m) 



The microscopic field equations 



c J- ^ 



AB 



by: Gab = 
tion n^ip = dV/dip - e^'^^'^'^f where, 




g5\^mix,^X,gAB)- 

5]) are now replaced 
, and the KG equa- 



p(m) 
-AB 



Sg 



AB 



AB ' 



K5 e -t-KgJ^j, Lj^jy 



gyiv 



i?-K5(V0) ^2 g5«5/3(v+0)r(™) 



(^'2 - 2V{v)) , 



9^ + «;,5/3e5«='5(^+^)T5('") = ^ _ ^" + 4i/^',(19) 

dip 



where f^f = 



g^iy(V(/))2/2 is the stress-energy 



tensor for the (massless) scalar field's component 
(/)(x). Besides, the consistency condition T^^^ = T^^'' /2 



in 



g'^'^T^r^) leads to Tg^"^ 



Srj^VS- Here we do 
not care about a concrete realization of a matter stress- 
energy tensor which fulfills the above conditions and just 
assume it can be found. The trivial case which meets 
these conditions is the radiation matter field. 

Under the ansatz ^ the following 5D solution is found: 



(20) 



where Vb = Z^"^ /2k\, and po is a non-vanishing integra- 
tion constant since, under the chameleon coupling, the 
mirror symmetry is not a symmetry of the field equations 
any more. The above 5D solution is consistent with the 
4D microscopic EFE-s: 



□0 = _3k5/3 g5K5/3(^+0) „2 



(21) 



Even if the potential V in (PU]) is unstable, what really 
matters to the dynamics of (/^(x, y) is not V itself, but 
the chameleon potential [l^ 



K// = Vb-2/l2(<^-<^0)2 + 



r e 



(22) 



is the chameleonic matter stress-energy tensor with trace 

and we have taken into ac- 



(m) 



(m) 



5 ~ g ^MN — ^5 

count the relationship 5g^^' = e^'^^l^'^S^S^^ Sg"^^ . In 



terms of the stress-energy tensor T^'^'' the resulting mi- 
croscopic EFE and the KGE are: 



Gab = K5 



^ AB + 4 R 



AB 



(18) 



respectively, or, if substitute the metric (IT51) into (|18l) 



4 < 0) this 

which solves the non- 
■21 



where Lp — p{y) is given by (j20[) . and we have defined 

r = r(x) = -3e5«^.'3'^'ri"VlO. If r > (t/"^ 
potential is an extremum at p 

algebraic equation: SKs/^rexp (bK^/Sp^,) = ^^'^{ip^ — pQ 
If ip^ > ipo + 1/5k5/3 the extremum is a minimum at 
Notice that, if assume /3 is a non-negative constant, 
then only for the positive branch of the solution (f20| the 
minimum of V can be found. In what follows we shall 
consider only the positive branch of (|20)) . The value p^, 



picks a position in the extra-space = -\/ Kg/S/i^ ((y9» — 
(^50)5 which meaning will be clear below. In the figure 
[5] the form of the potential is shown for arbitrarily 
chosen values of the parameters Vq, fj,, and /?. Different 
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FIG. 3: The 4D graviton's probability density »I'o(y) (solid 
line), the probability density of the massive states '^^{y) 
(dashed line), and the shifted probability density distribu- 
tion ^' shift (y) (doted line), are shown. The following values of 
the parameters have been chosen: (A, m, /3) = (0.5, 0.01, 0.5). 

constant values have been also assigned to the function r. 
However, since this is a function of the spacetime point, 
the 4D dynamics affects the existence and position of the 
minimum of the potential and, hence, the value y*. 

B. Mass Hierarchy 

The chameleon coupling shifts the maximum of the 
probability density distribution for the matter degrees 
of freedom with respect to the gravitational ones. To 
see this let us to average the first term in the RHS of 
Einstein's equations in (PT|) . 

^2^5/5«.(,+.o+lM) r Tl^^^^l^^^^y, (23) 

J — oo 

where ^>l^,,{y) = (y2A/A) e;2(y-?/o)VA^ ^ ^liv^V^): 
is the 'shifted' probability distribution function with j/q = 
5-\/3/3A/4. As it has been clearly illustrated in Fig|3l it 
is apparent that, while the probability density function 
for the gravitational degrees of freedom vE'o(j/) in Eq. pT]) 
is peaked at the origin j/ = 0, the one for the matter 
degrees of freedom ^g^ift peaked at y = j/Oi i- e., it is 
shifted in the extra-direction. 

Let us to do some juggling with the free parameters. 
First we write the integration constant lpq in terms of 
the values at the minimum of the effective chameleon po- 
tential ipo — — V^fiy^/n^ and substitute in Eg. ([23l) . 
Then we identify the position yo of the center of the 
shifted distribution function (^' shift) with the position 
y^, picked up by the minimum of V^ff: yo — y*. This 
choice forces the relationship (3 = (4/5\/3) 2/*/A. Since 
the integration constant (po is arbitrary, we can make yet 
another assumption on the constants: ip^, — a/3a'2/*/2k5 
(fiQ = —ifif, which greatly simplifies the resulting ex- 
pressions. After the above choice of the free parameters 
Eg.ip^ can be written in the following way: 



where "()^" means averaging over the shifted Gaussian 
distribution ^^(y) = '^ohj — V*)- The 4D macroscopic 
EFE-s - the ones which dictate the gravitational laws a 
4D observer like us measures - can be obtained through 
Gaussian averaging of the EFE-s (|2T|) . We obtain: 

n - 5K5P^ /rj.{Tn}\ , 2f{4>) 

Let us, for definiteness, to assume T^™' in the above 
equation can be written as r^™'(x, y) — (y) tIiT\'x.) 
=> rj'"' = rj™' (see the former section). Hence 




Accordingly, the averaged EFE-s 

G^. = 87rGAr(e5^^^^f^™)+f^^)), 

can be obtained if the following relationship between 
"Newton's gravitational constant" and the 5D gravita- 
tional coupling takes place: 

V z K5 

If assume that the fundamental 5D scale I/K5 ~ TeV, 
and take into account that the effective Planck scale 
Mpi = lO^^GcV, then, from the above equation one finds 

Hence, we do not need a large value of the constant pa- 
rameter /3 to explain the large hierarchy between the 
mass scales. Recall that, in harmony with string theory, 

VII. DISCUSSION AND CONCLUSION 

In the present paper we have addressed the issue 
of how to compute in an unambiguous way effective 
four-dimensional quantities within thick brane contexts 
[lOl [Tl| . Our proposal rests upon the following assump- 
tions: i) the background geometry is a warped one with 
a Gaussian warp factor and ii) the extra-coordinate 
y is regarded as a continuous randomly distributed pa- 
rameter. The resulting picture supports the following 
probabilistic interpretation. Fields which depend on the 
extra-coordinate are also randomly distributed (we call 
them as microscopic fields), their probability distribution 
density is a Gaussian with variance a = A/2, where A is 
the width of the brane. The corresponding 4D measur- 
able quantities coincide with their Gaussian averages. 

As in RS braneworlds 0, 0] the warped background 
geometry warrants that a bound 4D graviton exists and 
is capable of reproducing standard 4D laws of gravity in 
the brane. This time, however, due to the well-shaped 
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form of the analog quantum mechanics potential Uqm 
(see FigH]), even if the excited KK states form a con- 
tinuum, there exists an effective mass gap between the 
excited KK modes and the ground massless state. This 
is explained by the fact that the probability density of 
the KK states is a vanishing minimum at the brane lo- 
cation and, as we depart from the brane, it grows slower 
for the lighter KK modes. This means that the heavier 
KK excitations are closer to the brane than the lighter 
ones and, hence, can interact with the graviton with a 
greater probability. However, low energy states living in 
the thick brane need of an energy of the order of the mass 
of the KK modes to excite them. The interplay between 
probability of interaction and mass of the KK states is 
what generates the effective mass gap. 

Another ingredient of our approach is related with the 
way measurable 4D quantities are obtained. Here the ef- 
fective picture is generated by averaging the microscopic 
equations and quantities over the Gaussian distribution 
for the random parameter y to get the corresponding 
macroscopic equations and quantities, i. e., the ones a 
4D observer living in the thick brane is capable of mea- 
suring. We want to underline that the RS solution Q 
can not be recovered from our model under any circum- 
stances, i. e., the latter does not represent a thick brane 
generalization of the RS model as it is customary 
It is just a different phenomenological scenario. The fact 
that the solution (|6l7p supports a smooth thick brane 
configuration is a direct consequence of the (linear) grav- 
itational content of our model and of our choice of the 
way observable quantities are defined. 

In order to be able to explain the large hierarchy be- 
tween the TeV and Plack scales in our approach one 
needs to assume a non-minimal (chameleonic) interac- 
tion between the matter degrees of freedom living in 
the brane and the 5D scalar field. An interesting pic- 
ture arises. The resulting probability distribution func- 
tion for matter, which is shaped by is shifted in the 
extra-space with respect to the one for gravity, which 
is shaped by ^'q: "^l = ^'§(?/ — 2/»), where the position 

of (the center of) the matter probability distribution 
is determined by the minimum of the effective potential 
for (^(x,2/) = if ly) + (/)(x): K// ^ Vq - 2ii'^{ip - ipo)"^ + 
^g5K,5l3ip^ whenever the minimum exists. Besides, since 



T = -Se^'^^'^'^T^"''' /lO = t(x), the 4D spacetime dynam- 
ics of matter affects the 5D dynamics through modify- 
ing the position of the matter probability distribu- 
tion function. If assume, for instance, 4D matter in the 
form of a perfect fluid with energy density p(x) and pres- 
sure p = (7 — l)p (7 is the barotropic parameter), then 
t(x) cx (4 — 37) e^'^s^'^^'') p(x). For larger values of p the 
position y^, of the matter probability distribution func- 
tion 5'^ - the one which shapes the thick brane where the 
matter is trapped - is closer to y = where the graviton's 
probability distribution function is a maximum, while for 
smaller p-s it is farther away (see Figl2]). For a relativistic 
fluid with 7 = 2 (stiff- matter) , since t < 0, the minimum 
does not exist and the position of the matter probability 
distribution is not stabilized at all. Within a cosmological 
setting this would mean that, close to the big bang, when 
relativistic matter is expected to dominate the cosmic dy- 
namics, the position in the extra-space of the thick brane 
where matter lives (the one shaped by 5*^) is not stabi- 
lized. As long as non-relativistic matter starts dominat- 
ing, a minimum of Ve// is found and the matter brane's 
position is stabilized at some y*. As the expansion fur- 
ther proceeds and the matter energy density dilutes, the 
position of the matter brane gets farther and farther 
away from the brane at y = where the massless gravi- 
ton lives, which is shaped by ^q. This entails that the 
strength of gravitational interactions felt by matter con- 
fined to the thick brane at y^ weakens with the cosmic 
expansion. This happens after the radiation domination 
epoch (r = 0) where the scalar field's potential drives un- 
stable 5D dynamics. Before that stage, as said, perhaps 
relativistic matter dominated so that stabilization of the 
matter's thick brane position was not achieved neither. 
A more careful study, including detailed investigation of 
the 4D cosmic dynamics within the present approach, 
will lead to new interesting features not explored in the 
present paper. 
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